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Abstract 

We calculate the leading C(af ) contributions to the invariant mass distribution of top- 
quark pairs produced at the Tevatron and LHC, in the limit where the invariant mass 
of the ti pair approaches the partonic center-of-mass energy. Our results determine 
at NNLO in a s the coefficients of all singular plus distributions and scale-dependent 
logarithms in the differential partonic cross sections for qq,gg —> ti + X. A numerical 
analysis showing the effects of the NNLO corrections on the central values and scale 
dependence of the invariant mass distribution is performed. The NNLO corrections are 
found to significantly enhance the cross section and reduce the perturbative uncertainties 
compared to the NLO calculation. 



1 Introduction 



The top quark is the heaviest known particle in the Standard Model (SM) of fundamental 
interactions. Because of its large mass, it is expected to couple strongly with the fields re- 
sponsible for electroweak symmetry breaking, and the detailed study of top-quark properties 
is likely to play a key role in elucidating the origin of particle masses. In addition to observ- 
ables such as the top-quark mass and total inclusive top-quark pair production cross section, 
differential cross sections are also of interest. For instance, the ti invariant mass distribution 
can be used to measure m t [I], and in searches for physics beyond the SM. The presence of 
bumps in the smoothly decreasing ti invariant mass distribution would be a clear signal of an 
s-channel heavy resonance [2], which is predicted in many beyond-the-SM scenarios [3]. 

To date, thousands of top-quark events have been observed by two different experiments 
at the Fermilab Tevatron, and the top-quark mass has been extracted at the percent level. 
Searches for narrow width particles decaying into top-quark pairs have been pursued [Jj, and 
results for the top-pair invariant mass distribution were recently obtained from data collected 
by the CDF collaboration [5J. The experiments at the LHC are expected to observe millions 
of top-quark events per year already in the initial low-luminosity phase, bringing the study 
of top-quark properties into the realm of precision physics. In particular, the total inclusive 
top-quark pair production cross section is expected to be measured with a relative error of 5% 
to 10% at the LHC To make optimal use of the data requires equally precise theoretical 
predictions for the measured observables. 

Theoretical calculations in QCD rely on the factorization formula for the differential cross 
section, which is of the form 



where the symbol <g> stands for a convolution. The hard-scattering kernels Cy are related to 
the partonic cross section and can be calculated as a series in a s , whereas the /j are parton 
distribution functions (PDFs) for the partons % = q, q, g in the incoming hadrons and must be 
extracted from data. Current theoretical predictions are based on next-to-leading order (NLO) 
calculations of the total cross section [7HTU] and differential distributions [TT|[T2] in fixed-order 
perturbation theory, supplemented with soft gluon resummation to next-to-leading-logarithmic 
(NLL) order [T3HTT] . The NLO computations suffer from theory uncertainties larger than 
10%, both for Tevatron and LHC center-of-mass energies. These uncertainties are due to our 
imperfect knowledge of the parton distribution functions, and also to the truncation of the 
perturbative series in the strong coupling constant, which introduces a dependence on the 
unphysical renormalization and factorization scales into physical predictions. This theoretical 
uncertainty is typically reduced by including more terms in the perturbative series, and for this 
reason the calculation of the differential partonic cross section to next-to-next-to-leading order 
(NNLO) has been an area of active research [T8H22] . However, due to the complexity of the 
diagrammatic calculations, complete NNLO results are not yet available. In the absence of full 
results, there has been much recent activity in NNLO calculations and next-to-next-to-leading- 
logarithmic (NNLL) resummation for the total cross section near production threshold 
but the analogous results for differential cross sections remain unknown. 

The goal of this Letter is to present a subset of the NNLO corrections to the hard-scattering 
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kernels for the ti invariant mass distribution at hadron colliders. In particular, we focus on 
the region where the invariant mass M of the top-quark pair approaches the partonic center- 
of-mass energy Vs, and compute NNLO corrections to the hard-scattering kernels for the qq 
and gg channels at leading order in 1 — z, where z = M 2 /s — > 1 in the threshold region. The 
leading-order term in this threshold expansion for the differential cross section is equivalent to 
the virtual-soft approximation, and is written in terms of singular plus distributions and delta 
functions in the variable 1 — z. Our results determine the coefficients of all plus distributions of 
the form [ln n (l —z)/(l — z)] + , as well as all /i-dependent pieces multiplying the 5(1 — z) term, 
whereas a remaining piece of the delta-function coefficient is left undetermined. The basis 
for our calculations is a factorization formula for the hard-scattering kernels in the partonic 
threshold region, and is described in Section |5J In this region, the hard-scattering kernels fac- 
torize into products of hard functions, related to virtual corrections, and soft functions, related 
to real emission in the soft limit [13]. These functions satisfy certain renormalization-group 
equations, which determine their dependence on the scale p. By calculating the hard and soft 
functions at one-loop order, and using results for the two-loop anomalous dimensions recently 
derived in [2"5H3T| (see also J32H2H]), the /i-dependent logarithms in the hard and soft functions 
can be determined exactly to NNLO using the renormalization group. The logarithms in the 
soft function are of the form ln(s(l — z) 2 /p 2 ), and uniquely determine the coefficients of the 
[ln n (l — z)/(l — z)] + distributions. This is similar in spirit to the calculations of [16] for the soft 
corrections to the NNLO differential cross section using threshold resummation techniques in 
Mellin space, but goes beyond those results by completely determining the coefficient of the 
[1/(1 — z)] + distribution, which is sensitive to process-dependent two-loop anomalous dimen- 
sions. In Section |3] we perform a short numerical analysis illustrating the impact of the NNLO 
corrections on the central values and scale dependence of the ti invariant mass distribution. 



» P = J ^h ( 2 ) 



2 NNLO corrections at threshold 

In this section we discuss our method for determining NNLO corrections to the differential 
PP (pp) tt + X cross section in the threshold region. To describe the threshold region we 
introduce the variables 

M 2 M 2 Am 

~M 2 

where M is the invariant mass of the ti pair, a/s is the hadronic center-of-mass energy, and 
\/I is the partonic center-of-mass energy. We define the threshold region as the limit z — > 1, 
with p = 4m 2 /M 2 a generic 0(1) variable, which is however not too close to unity. Note that 
this is different from the threshold limit p — > 1 often used for the total cross section. 

The calculation of the differential cross section in the threshold region is greatly simplified, 
since there is little phase space available for real gluon emission. Hard emissions are suppressed 
by powers of 1 — z, and the partonic scattering process is dominated by hard virtual corrections 
and the real emission of soft partons, which can be calculated in the eikonal approximation 
and exponentiate into Wilson lines. The phase space for these processes is effectively two- 
body, and the fully differential cross section can be written in terms of the kinematic variables 
entering the Born- level result. 
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At the Born level the differential cross section receives contributions from the quark- 
antiquark annihilation and gluon fusion channels 



q(pi) + q(p 2 ) -»■ t{p 3 ) + t(pi) , 

(3) 

g(pi) + g(p2) -> t(p 3 ) + t(p A ) . 

The partonic cross section is a function of the kinematic invariants 

s = (pi+P2) 2 , h = (pi - Pz) 2 - ml , Mi = (pi - p 4 ) 2 - m 2 t , (4) 

and momentum conservation implies s + 1\ + "Ui = 0. The hadronic cross section is obtained 
from the partonic one by convoluting with PDFs. The fully differential cross section involves 
three variables, which can be chosen, for instance, as the invariant mass and rapidity of the ti 
pair, and the scattering angle 8 between p\ and p 3 in the partonic center-of-mass frame. The 
kinematic invariants in (J3J) can be written in terms of these variables according to 

2 ]\/T2 



s= , h = -— (l-/9cos0), p=^l^~p. (5) 

Less differential results are derived by integrating the triply differential rate over the appro- 
priate phase space. 

In this Letter we will be interested in the invariant mass distribution of the ti pair. We 
shall write this distribution at threshold as 

da thvcs 8n/3 f 1 dz r /r \ /r \ ] 

~dW~ = 3iM J T r" w v Cg9 ^ z ' M ' m *' ^ + ffqq W ^) Cq ^ z ' M ' mt1 ^\ + °^ ~ ^ ' 

(6) 



where we have introduced the parton luminosity functions 

f x dx 

ffij(y^)= —fi(x,fj,)fj[y/x,fj,). (7) 

Jy X 

For simplicity, we set the factorization and renormalization scales equal and refer to them as 
jj,. The luminosities for qq are understood to be summed over all species of light quarks, and 
we have made explicit that mixed channels such as qg are power suppressed by 1 — z in the 
threshold region [13j. The hard-scattering kernels Cy factorize into hard functions, related to 
virtual corrections, and soft functions, related to soft real emissions [13]. The result is 

Ci(z, M,m t , /i) = J <icos#Tr H^M, m t , cos 9, /i) Si(y/l(l — z), M, m t , cos 8, fx) . (8) 

We have used a label % = qq, gg to distinguish the different partonic channels, but will suppress 
it in the formulas below, referring instead to a single coefficient function. Note that the hard 
and soft functions are matrices in a color decomposition of the QCD amplitudes. We have 
used s = M 2 /z ~ M 2 everywhere except in the first argument of the soft function, since we 
expect that the exact results contain logarithms of the form m(s(l — z) 2 / /i 2 ), as in the case of 
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Drell-Yan production [27]. In J3S], we will give details on the derivation of the factorization 
formula (jHJ) in soft-collinear effective theory. 

The perturbative expansion of the hard-scattering kernel can be written as 



C = ai 



4:71 \AtxJ 



(9) 



The purpose of this Letter is to present new results for the NNLO coefficient C (2) . As men- 
tioned in the introduction, our derivation relies on the renormalization-group equations for 
the hard and soft functions. The soft function contains singular plus distributions, and its 
renormalization-group equation is non-local. It is convenient to avoid this complication by 
working instead with the Laplace transform of the soft function [39JSQ] , defined as 



POD 

s(L, M, m tl cos 9, fi) = / duo e" aoJ S(oj, M, m t , cos6*, fi) , a 

Jo 



(10) 



Whereas the soft function contains distributions, the Laplace-transformed function is a polyno- 
mial in its first argument and satisfies a local evolution equation. To rewrite the hard-scattering 
kernel in terms of this function, we introduce the notation 



c(d v , M, mt, cos 9, /i) = Tr H(M, m t , cos 9, /i) M, m t , cos 9, /i) 



where d v is a differential operator with respect to an auxiliary variable rj. The factorization 
formula for the hard-scattering kernel then takes the form [38] 



C(z,M, m t ,fj) 



d cos 9 c(d v , M, m t , cos 9, n) 



(M\ 2v e~ 2 ^ 



2,i 



(12) 



»?=o 



To evaluate the above formula in terms of distributions in the variable z, one must use analytic 
continuation to regulate the divergence at z — > 1, take derivatives with respect to the auxiliary 
variable i], and then set rj — > 0. This procedure yields results in the form of plus distributions 
and delta functions given in (T2~T|) below. 

From the above discussion, we conclude that it is sufficient to focus on the calculation of c 
at NNLO. Defining its perturbative expansion in analogy with (Q, and using the expansions 



H 



H w + ^ H (i) + (^.\ 2 H (2) + 

47T V47T/ 
47T V47T/ 



(13) 



for the hard and soft functions (El), the result at NNLO reads 



^{d v ,M,m t , cos 9,n) = Tr [H ^ I<°>] + Tr [H ® I< 2 >] + Tr [ffW^] 

4 

= Y,c?\M,m t ,cos9,n)di l . 

3=0 



(14) 



4 



In the second line we have defined expansion coefficients multiplying explicit powers of d rj . We 
see that obtaining a complete result for the threshold expansion at NNLO would amount to 
calculating the hard and soft matrices to this order. Although the NLO coefficient c^ 1 ^ can be 
extracted from known results for the virtual corrections and real emissions in the soft limit, 
the hard and soft matrices themselves are so far unknown, except for in the absolute threshold 
limit p — > 1 A main result of our work is the calculation of these functions at NLO, for 

both the qq and gg channels; details of the calculation will appear in [38]. As for the NNLO 
corrections, we will show in what follows that it is possible to extract all of the coefficients 

(2) 

except for Cq by using the renormalizat ion-group equations for the hard and soft functions, 
along with the one-loop results and the two- loop anomalous dimensions from [29"H3"T] . The 
renormalizat ion group also determines the /^-dependent logarithms in this remaining term, 
but to obtain the scale-independent piece would require to evaluate at NNLO the virtual 
corrections and soft real emissions. 

We now sketch how to calculate the coefficient in (|T^|) using the evolution equations for 
the hard and soft functions. The hard function satisfies 



— — H(M, m t , cos 9,fi) = T H H(M, m t , cos 9, (i) + H(M, m u cos 9, u) T H , (15) 

CL 111 ^Jj 

where Th is given by the matrices T q q or T gg from (31] . The anomalous dimensions to two- loop 
order can be decomposed as 

M 2 ■ V .hi 



T H = T cusp (a s ) ( In — - iir ) 1 + 7 (M, m t , cos 9, a s ) , (16) 



where the object 7 is defined through a comparison with the explicit results of [31J, and 
Tcusp is equal to CV7cus P for qq and CaJ C us P for gg, with 7 cusp the universal cusp anomalous 
dimension. The evolution equation for the soft function follows from that for the hard functions 
and PDFs, along with RG-invariance of the cross section. The result is 



s In — — , M, mt, cos 9, jj, 



dlia.fi \ \i 



T\ s ( In — — , M, m t , cos 9, a) + s I In — — M, m t , cos 9, a) T s 
\ fi 2 J \ fi 2 J 



(17) 



where the soft anomalous dimension is given by 

/ M 2 \ 
T, = - r cusp (a s ) fin— -in J + 2^(a s 



1 - j h (M, m t , cos 9, a s ) 



The object 7^ is defined through the large- a; limit of the Altarelli-Parisi splitting functions, 
which reads 

P(x) = g^gP^f) + 2 /(« s ) 6(1 - X ) . (19) 
(1 - x)+ 

Since the anomalous dimensions for the hard function and PDFs are known to two-loop order, 
we can calculate that of the soft function to this order. 
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Given the evolution equations, the anomalous dimensions, and the one-loop hard and soft 
matrices, it is a simple matter to deduce a general expansion of the form 



2n 



M 2 

Tr \H in) s (0) l (M, m t , cos 9, a) = V h^(M, m t , cos 9) ln J — , 



2n 



(20) 



Tr [H {0) s {n) ] (d v ,M,mt, cos 9) = ^ {M,m t , cos 9) d J v 



and to calculate all of the coefficients up to and which are left undetermined. Along 
with the Tt[H^s^] contributions, these determine c{\... } cf^ in (JSJ), as well as the /in- 
dependent piece of c[, 2 \ 

The hard-scattering kernel C is obtained by evaluating (112)) . In doing so, powers of in 
c are converted into delta functions and plus distributions in the variable z. The final result 
can be written as 



C {2 \z,M,m t ,fx) 



d cos 9 I Di 



+ D 1 



ln(l - z) 
1 - z 



ln d (l - z) 
l-z 



Dr 



+ D 2 



ln 2 (l - z) 
1 - ^ 



(21) 



1 - z 



+ C 5{l-z) + R{z) }. 



Explicit results for the coefficients in both the gg and qq channels as functions of the variables 
(M, mt, cos#, /i) are given in the computer program attached with the electronic version of 
this article. We stress that the coefficients D Q , . . . , _D 3 are completely determined from our 
calculations. The coefficient Cq, on the other hand, is correct only in its /i-dependence. It also 
receives scale-independent contributions related to the translation from c to C, and from the 
product of NLO hard and soft functions, which we discard for consistency. To determine the 
remaining scale-independent piece would require a complete NNLO calculation of the hard 
and soft matrices. Note that, because this piece is unknown, the expansion of Cq in terms 
of ln(M 2 //i 2 ) considered so far is not unique: expansions in terms of ln(m 2 //i 2 ) (or other 
choices) still satisfy the same evolution equations. In our phenomenological study in the next 
section we shall consider these two different expansions as one indication of systematic errors 
related to the missing, scale-independent piece of Cq. The function R(z) is finite in the limit 
z — > 1 and is not determined from the leading-order threshold expansion. However, it receives 
a contribution related to the z~ v factor in ( fl2|) . We choose to keep these contributions when 
evaluating the NLO corrections, but drop them from the NNLO ones. Our results for Di, 
D 2 , D 3 , and the /i-dependent terms agree with [IB], while the exact result for D contains 
process-dependent two-loop effects and is new. 

The coefficients Di and Cq do not depend on the variable z. It is therefore possible to 
convolute the plus distributions and delta function with the luminosities as in (J6j). For the 
phenomenologically interesting values of 0.03 < r < 0.3 discussed in the next section, we can 
then make a few general comments concerning the hierarchy of the different plus-distribution 
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and delta-function terms. The contribution of the leading [ln 3 (l— z) / (1— z)} + distribution is the 
largest in all cases. The [ln 2 (l — z)/(l — z)]+ distribution is the second largest, its contribution 
being about 1/3 that of the leading one. The contribution of the delta-function is about 10- 
20% of the leading plus distribution, and is of comparable size to the [ln(l — z)/(l — z)] + and 
[1/(1 — z)] + distributions. 

3 Phenomenological applications 

In this section we study the phenomenological implications of our results. The main goal is 
to show the impact of the NNLO corrections on the central values and perturbative uncer- 
tainties of the ti invariant mass distributions at the LHC and Tevatron. However, given that 
our NNLO results are incomplete, it is important to clarify under what conditions they are 
expected to give a good approximation to the full results. While it is difficult to assign a 
systematic uncertainty to an unknown correction, we can examine in parallel how well the 
same approximation works at NLO, and use this as an indication of its expected validity at 
NNLO. 

At NLO we can compare the following three calculations of the invariant mass distribution: 

1. The full result, given by the Monte Carlo program MCFM jUJ. 

2. The full threshold expansion, which is equivalent to the full result at leading order in 

1 -2. 

3. The approximate threshold expansion, which is equivalent to the full result at leading 
order in 1 — z for the plus distributions and //-dependent terms, but not in the /x- 
independent piece of the delta-function term. Here we also have the choice of writing 
the logarithms in the delta-function coefficients in terms of ln(m 2 //i 2 ) or ln(M 2 /// 2 ). 

4. The approximate threshold expansion, where the delta-function term is neglected en- 
tirely. 

To the extent that the second, third, and fourth options agree well with the full result at NLO, 
it would seem plausible that our NNLO result, which is limited to the third and fourth options, 
might also agree well with the full NNLO result, which is unknown. We thus investigate these 
different approximations to the invariant mass distributions for a range of M and /i, and study 
under what circumstances they show a reasonable agreement. The naive expectation would 
be that they agree only at very high values of the invariant mass, since then r = M 2 /s — > 1 
and the integrand in is needed only in the z — > 1 limit. Moreover, the integrals over 
the plus distributions in ( 12 ip are more singular at r — > 1 than the unknown portion of the 
correction proportional to the delta function. However, due to the rapid fall-off of the PDFs 
at large x, the distribution at very high r is difficult to measure. The most interesting region 
for phenomenology would be from M ~ 2m t to around 1 TeV at the Tevatron and up to 
several TeV at the LHC; this region of the invariant mass corresponds to r < 0.3. In order 
for our results to be a good approximation in this region, it is necessary that the luminosity 
functions ffij(r/z, /i) fall off so fast for rjz — > 1 that only the largest values of z give significant 
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M (GeV) M (GeV) 

Figure 1: Ratio K = (da / dM) / (da LO (p = 2m t )/dM) versus the tt invariant mass M for the 
Tevatron (left) and the LHC with yfs = 10 TeV (right). The binned histograms show the exact 
NLO results from the Monte Carlo program MCFM [UJ, while the solid curves are obtained 
from the threshold expansion at NLO. In each case the upper lines correspond to /i = mt and 
the lower ones to /i = 4mt. 

contributions to the integrand in (JBJ). In that case, a dynamical enhancement of the partonic 
threshold region may occur, even if r is not close to unity [3~7 l l4"2 | l4*3]. The agreement of 
the three different approximations at moderate values of r is thus a measure of whether this 
dynamical enhancement actually occurs. 

With these points in mind, we begin by comparing the exact result and the complete 
threshold expansion at NLO. The results are shown in Figure [TJ Here and throughout the 
analysis we usejhe MSTW2008 NNLO PDFs 03], take a s (M z ) = 0.117 with three-loop 
running in the MS scheme with five active flavors, and use m t = 173.1 GeV in the pole 
scheme. We use the same set of PDFs at each order of perturbation theory so as to better 
illustrate the size of the different perturbative contributions to the hard-scattering kernels CV,-. 
We choose \x = 2m t as the central value and vary it between m t and 4m t . For the case of the 
Tevatron, we see that the threshold expansion is in good agreement with the full result even 
at values M ~ 2mt, where r w 0.03, and gets progressively better at higher values of M. The 
situation is only slightly worse for the LHC with y/s = 10 TeV, even though the value of r 
for a given M is about 25 times smaller. We also note that at LHC energies the threshold 
expansion matches the exact result better at the higher value of \i = 4m t than at fi = m t . In 
what follows, we shall also explore the scale choice \x ~ M. Since such a choice is not practical 
in the Monte Carlo program MCFM, from now on we shall use the exact threshold expansion 
to represent the full NLO result, keeping in mind the good agreement seen in Figured! 

We next investigate in more detail the approximate threshold expansions (options 3 and 
4 above) at NLO and NNLO. In Figures [2] and [3] we compare the full threshold expansion 
at NLO with the different approximate expansions at the same order, and also show the 
approximate expansions at NNLO (these are generated by adding the NNLO corrections to 
the full threshold expansion at NLO). In Figure [2] we show the differential cross section as 
a function of \x for two different values of M. We notice that the different approximations 
at NLO show better agreement with the full threshold expansion for higher values fi ~ M, 
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Figure 2: Scale dependence of the tt invariant mass distribution at the Tevatron (upper plots) 
and LHC (lower plots), for two different values of M. The dotted, dashed, and solid lines 
show results obtained using the threshold expansion at LO, NLO, and NNLO, respectively. 
The long-dashed black lines are obtained using the full threshold expansion at NLO. The 
light-grey (dark-grey) lines correspond to the approximate threshold expansion (option 3), 
in which the logarithms in the delta-function coefficients are written in terms of ln(m 2 //z 2 ) 
(ln(M 2 //x 2 )), and the black lines to the threshold expansion where the delta-function term is 
dropped altogether (option 4). 

especially at LHC energies, where the NLO approximations at lower values of fi differ greatly 
from the full results. At NNLO, the different approximations in options 3 and 4 do not differ 
substantially from one another, when one considers their highest and lowest values in the 
range M/2 < // < 2M, which is how perturbative uncertainties are typically estimated. In 
Figure |3] we show results for the same approximations at NLO and NNLO, but this time as a 
function of M. Given the better agreement at higher /i observed in the previous two figures, 
we have made the default choice \i = M. We observe that at NLO the approximate threshold 
expansions recover a significant portion of the exact NLO correction, both at the Tevatron 
and LHC. It is not unreasonable to expect that the same is true of our approximate NNLO 
results, although they are clearly no substitute for a complete NNLO computation. 

In Figure H] we show the invariant mass spectrum as a function of M, with bands repre- 
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Figure 3: Ratio K = (da/dM)/(da LO (fi = M)/dM) versus the ti invariant mass M for the 
Tevatron (left) and LHC (right), for the default scale choice \i = M. The dashed and solid 
lines show results obtained using the threshold expansion at NLO and NNLO, respectively. 
The black-dashed lines are obtained using the full threshold expansion at NLO. The light-grey 
(dark-grey) lines correspond to the approximate threshold expansion (option 3), in which the 
logarithms in the delta-function coefficients are written in terms of ln(m 2 //i 2 ) (ln(M 2 //z 2 )). 



senting the uncertainty associated with scale variations in the range M/2 < fi < 2M. For 
the NNLO result, we have displayed only the approximate threshold expansion in which the 
logarithms in the delta-function coefficient are written in terms of ln(M 2 //i 2 ); those in the 
other two NNLO approximations are similar, and very nearly contained within this band. At 
both the Tevatron and the LHC we note a reduction of scale uncertainty upon including the 
NNLO corrections, which is the expected benefit of including more terms in the perturbative 
expansion. However, we must also observe from Figure H] that at NLO the perturbative uncer- 
tainty estimated by varying the scale in the range M < \i < 2M in the approximate threshold 
expansions is actually smaller than in the complete threshold expansion at both M = 400 GeV 
and M = 1000 GeV at the LHC. We thus cannot rule out that the very small uncertainty of 
the NNLO result at small values of M is also an underestimate. 



4 Conclusions 

We calculated 0(af) contributions to the it invariant mass distribution at hadron colliders. 
The calculation was based on an effective field-theory approach to the factorization of the per- 
turbative hard-scattering kernels into matrix-valued hard and soft functions in the partonic 
threshold region s > M 2 . At the technical level, it involved calculating the hard and soft func- 
tions at NLO, and then determining at NNLO the logarithmic /i-dependence of these functions 
using the renormalization group and recent results for two-loop anomalous dimensions. Our 
computations yield exact results for the coefficients of all singular plus distributions and in- 
dependent logarithms in the differential partonic cross section at NNLO, in the limit where the 
invariant mass of the ti pair approaches the partonic center-of-mass energy. They agree with 
those previously obtained in [16] using a Mellin-space approach to soft gluon resummation, 
but go beyond them by uniquely determining the coefficient of the [1/(1 — z)] + distribution, 
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Figure 4: Top: Ratio K = (da / dM) J (do (fi = M)/dM) versus the ti invariant mass for the 
Tevatron (left) and LHC (right), with the factorization scale varied in the range M/2 < \l < 
2M. The light bands between the dotted lines correspond to the threshold expansion at LO, 
the medium bands between the dashed lines to NLO, and the dark bands between the solid 
lines to NNLO. Bottom: Same as above, but for the invariant mass distribution dajdM. 



which contains process-dependent two-loop effects related to anomalous dimensions of the soft 
and hard functions. To obtain the missing scale-independent piece of the delta-function term 
would require the exact two-loop calculation of these functions. 

In Section [3] we performed a numerical study of the ti invariant mass distribution at Teva- 
tron and LHC energies. In both cases, we noted the expected decrease in scale uncertainty 
after including the NNLO corrections, and also noticeable changes in the central values, espe- 
cially at higher values of the invariant mass. Given that our NNLO result is still incomplete, 
even in the threshold region, part of our analysis was focused on determining how well the 
analogous approximation works at NLO. We observed that at NLO the full result in QCD 
and its exact threshold expansion agree fairly well. The same was true of the approximate 
threshold expansion at NLO, where the scale-independent pieces of the delta-function term are 
neglected. This provides some evidence that the steep fall-off of the PDFs at large x induces 
a dynamical enhancement of the partonic threshold region, which would mean that the terms 
we have calculated at NNLO provide a useful approximation to the full result. 

Our finding that the NNLO perturbative corrections, which we have calculated in this 



11 



Letter, have quite an important impact in the region of large invariant masses, as seen in 
Figure HJ provides a strong motivation for improving this estimate with more sophisticated 
calculations. We expect that the scale dependence indicated by the width of the various bands 
in the figure can be reduced significantly by resumming the partonic threshold logarithms to 
all orders in perturbation theory. Based on the results presented here and in [SUET], such a 
resummation is now possible at NNLL order. It will be discussed in a forthcoming paper |38j . 
Ultimately, however, only a complete NNLO calculation of the spectrum will provide us with 
a complete picture of the true residual uncertainty in the theoretical predictions. 
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